Presentation
About things

# Define a server for the Shiny app
function(input, output) {
  
  # Fill in the spot we created for a plot
  output$graph <- renderPlot({
    # Render a barplot
  })
}

import numpy as np
import matplotlib.pyplot as plt

r = np.arange(0, 2, 0.01)
theta = 2 * np.pi * r
fig, ax = plt.subplots(subplot_kw={'projection': 'polar'})
ax.plot(theta, r)
ax.set_rticks([0.5, 1, 1.5, 2])
ax.grid(True)
plt.show()

library(dplyr)
library(ggplot2)

# ice cream sales data
sales <- c(5,5,10,5,15,20,35,40,35,20,10,5)
type <- c("not summer", "not summer", "not summer", "not summer", "summer", "summer", "summer", "summer", "summer", "not summer", "not summer", "not summer")
months <- c("Jan-23","Feb-23","Mar-23","Apr-23","May-23","Jun-23","Jul-23","Aug-23","Sep-23","Oct-23","Nov-23","Dec-23")

# make df
df_ices <- data.frame(sales=sales, months=months, type=type)

summer_months <- c("May-23","Jun-23","Jul-23","Aug-23", "Sep-23")

# get sum of sales for summer
sum_sales <- df_ices %>%
  filter(months %in% summer_months)

# get sum of sales for not summer
not_sum_sales <- df_ices %>%
  filter(!(months %in% summer_months))

# compare numbers and get p-value
sum_sales <- sum(sum_sales$sales)
not_sum_sales <- sum(not_sum_sales$sales)


# Hypothesis to test: sales are higher in the summer

proa<-prop.test(not_sum_sales,sum_sales, p=0.5, correct = FALSE)
#proa
proa <- list(proa)
proa <- as.data.frame(matrix(unlist(proa),nrow=length(proa),byrow=TRUE))
colnames(proa) <- c("X-squared","df","p-value","prop 1","prop 2","95pct conf int1","95pct conf int2","hypothesis","test name","vectors")

########## graph the p-value ##########################
graph <- rbind(proa$`p-value`)
graph <- as.data.frame(as.numeric(graph))
colnames(graph) <- c("p-value")
graph$xaxis <- c("Summer compared to not summer")

display<- ggplot(graph,aes(x=xaxis,y=`p-value`))+
  geom_bar(aes(fill=xaxis),position="dodge",stat="identity")+
  geom_text(aes(label=round(`p-value`,2)),vjust=-1.0,size=3)+
  scale_fill_manual(values=alpha(c("#CD2456","#14022E","#CD1234")))+
  #ylim(0,max(graph$`p-value`)+1)+
  theme_classic(base_size=10)+
  theme(legend.position = "none")+
  labs(x="AB Test",y="P-value",title="P-value for AB test")